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Abstract. In this article, we study an elastic manifold in quenched disorder in the limit of 
zero temperature. Naively it is equivalent to a free theory with elasticity in Fourier-space 
proportional to k 4 instead of k 2 , i.e. a model without disorder in two space-dimensions 
less. This phenomenon, called dimensional reduction, is most elegantly obtained using 
supersymmetry. However, scaling arguments suggest, and functional renormalization shows 
that dimensional reduction breaks down beyond the Larkin length. Thus one equivalently 
expects a break-down of supersymmetry. Using methods of functional renormalization, 
we show how supersymmetry is broken. We also discuss the relation to replica-symmetry 
breaking, and how our formulation can be put into work to lift apparent ambiguities in standard 
functional renormalization group calculations. 
Dedicated to Lothar Schafer at the occasion of his 60th birthday. 



1. Introduction 

The statistical mechanics of even well-understood physical systems subjected to quenched 
disorder still poses major challenges. For a large class of these systems, as e.g. random- 
field models or elastic manifolds in quenched disorder, an apparent simplification appears: 
Supposing that all moments of the disorder are finite, one can show that all correlation 
functions in the disordered model, in the limit of zero temperature, are equivalent to those 
of the pure system at finite temperature in two space-dimensions less, at a temperature 
proportional to the second moment of the quenched disorder. This phenomenon is called 
dimensional reduction (DR) [ 1 1. The most elegant way to prove it is to use the supersymmetry 
approach Q, as we will detail below. However, one also knows that dimensional reduction 
gives the wrong result at large scales, more precisely at scales larger than the Larkin length. 
The latter is obtained from an Imry-Ma type argument due to Larkin, balancing elastic energy 
and disorder, as we detail below. For a d-dimensional elastic manifold in quenched disorder, 
the elastic and disorder energy are 

E A [u] = J & d x l -(yu(x)f , EboN = J d d xV(x,u(x)) . (1) 
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For d = 1, these are polymers, for which a lot is known fl3); for d = 2 membranes; and 
for d = 3 elastic crystals, as e.g. charge-density waves. For simplicity we consider disorder 
which at the microscopic scale is Gaussian and short-ranged with second moment 



V(x, u)V{x', u') = 8 d (x - x')R{u - v!) . (2) 

Long-range correlated disorder R(u) is possible, and leads in general to a different 
universality class. This will play no role in the following. The most important observable 
is the roughness exponent (, which describes the scaling of the 2-point function 

[u(x) — ^(a;')] 2 ~ \x — x'\ 2C - . (3) 

The Larkin argument compares, as a function of system size L, elastic energy E c \ ~ L d ~ 2 
and disorder energy E DO ~ L d l 2 to conclude that in dimensions smaller than four, disorder 
always wins at large scales, leading to an RG-flow to strong coupling (in a way to be made 
more precise below). This suggests that the dimensional reduction result, derived below via 
the Supersymmetry method, 

-R"(0) , 4 d 

u kU- k = Cdr=^— (4) 

will become incorrect below four dimensions. 

2. The functional RG treatment 

In this section we review some important points of the functional RG treatment, which 
will facilitate the derivation of the corresponding formulas in the supersymmetric treatment. 
Functional RG was first introduced in [@J and pioneered for the problem at hand in 
13 □ El 0, to cite the earliest contributions. Important improvements ifTUl fTTl fHl fT3l flU 
GlEBGIlElGlIliEni^ have been obtained by several 

authors, see lf3TTl for a more detailed introduction and review. 

The Larkin argument suggests that four is the upper critical dimension and that an e- 
expansion Q with 

e = 4 - d (5) 

about dimension four is possible. Taking the dimensional reduction result © in d = 4 
dimensions tells us that the field u is dimensionless. Thus, the width a = —R"(0) of the 
disorder is not the only relevant coupling at small e, but any function of u has the same 
scaling dimension in the limit of e = 0, and might thus equivalently contribute. The natural 
consequence is that one has to follow the full function R(u) under renormalization, instead 
of just its second moment R"(0). Such an RG-treatment is most easily implemented in the 
replica approach: The n times replicated partition function becomes after averaging over 
disorder 



exp 



a=l a=l 

(1 n 1 11 f \ 

--^£ e iK] + — j d d xR(u a {x)-u b (xyj\ . (6) 

a=l o,6=l ^ / 
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Perturbation theory is constructed along the following lines (see ifTTl |2"61 for more details.) 
The bare correlation function, graphically depicted as a solid line, is with momentum k 
flowing through and replicas a and b 



T5, 



ab 



k 2 



(7) 



The disorder vertex is 



X 



* b 



J2R(ua(x)~U b (x)) . (8) 



a.b 



The rules of the game are to find all contributions which correct R, and which survive in the 
limit of T — > 0. At leading order, i.e. order R 2 , counting of factors of T shows that only 
the terms with one or two correlators contribute. On the other hand, ^ a b R(u a — u b ) has two 
independent sums over replicas. Thus at order R 2 four independent sums over replicas appear, 
and in order to reduce them to two, one needs at least two correlators (each contributing a 5 ab ). 
Thus, at leading order, only diagrams with two propagators survive. These are the following 
(noting C{x — y) the Fourier transform of 1/ k 2 ): 



v 



b * * b 



* b 



R"(u a (x) - u b (x))R"(u a (y) - u b (y))C(x - y) 2 (9) 



R"(u a (x) - u a (x))R"(u a (y) - u b (y))C(x - y) 2 . (10) 



In a renormalization program, we are looking for the divergences of these diagrams. These 
divergences are localized at x = y, which allows to approximate R"(u a (y) — u b (y)) by 
R"(u a (x) - u b (x)). The integral J x _ y C(x - y) 2 = f k {k 2^ m 2 )2 = ^ (using the most 
convenient normalization for J k ), is the standard 1-loop diagram, which we have chosen to 
regulate in the infrared by a mass, i.e. physically by a harmonic well which is seen by the 
manifold. 

Note that the following diagram also contains two correlators (correct counting in powers 
of temperature), but is not a 2-replica but a 3-replica sum: 




(11) 

* * c 

Taking into account the combinatorial factors, and a rescaling of R (which remember 

has dimension e for a dimensionless field u) as well as of the field u (its dimension being the 

roughness exponent Q, we arrive at 

d 1 
- m—R(u) = (e - 4C)R(u) + (uR'(u) + -R"(u) 2 - R"(u)R"(0) . (12) 
dm 2 

Note that the elasticity does not get renormalized due to the statistical tilt symmetry u(x) — > 
u(x) + ax. 

The crucial observation is that when starting with smooth microscopic disorder, 
integration of the RG-equation leads to a cusp in the second derivative of the renormalized 
disorder at the Larkin-length, as depicted on figure [T] This can easily be seen from the flow- 
equation of the fourth derivative (supposing analyticity), which from (fT2l is obtained as 

_ m JL#'"( ) = € R""(0) + 3R""(0) 2 . (13) 
dm 
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-R"(u) A 




-R"(u) A 



renormalization 

59- 




Figure 1. Change of —R"(u) under renormalization and formation of the cusp. 



(Note that this explains also the appearance of the combination e — 4£ in (fT^t ). This equation 
has a singularity R""(0) = oo after a finite renormalization time, equivalent to the appearance 
of the cusp, as depicted on figure [T] After that dimensional reduction © is no longer valid. 
This can most easily be seen from the flow of -R"(0): Deriving (fl2T) twice w.r.t. u, and then 
taking the limit of u — > leads to 

- J?-R"(0) = (e - 2OR"(0) + R"'(0 + ) 2 . (14) 
dm 

In the analytic regime R"'(0 + ) = 0, such that the fixed-point condition — -£^R"(0) = 
implies ( = | = after appearance of the cusp, -R'"(0 + ) ^ 0, thus ( has to change. 

This analysis can be continued to higher orders. Let us cite some key results at 2-loop 
order lfT4ll26l . for which the RG-equation reads 

- m^-R(u) = (e - 4C) R(u) + (uR'(u) + -R"(u) 2 - R"(u)R"(0) 

dm 2 

+ ~ (R"(u) - R"(0)) R"\u) 2 - ^R"'(0 + ) 2 R"(u) . (15) 

Different microscopic disorder leads to different RG fixed points. The latter are solutions of 
equation (fT5l) . with —m-£^R(u) = 0; it is important to note that given a microscopic disorder, 
the exponent solution of (fT3t is unique. For random-bond disorder (short-ranged potential- 
potential correlation function) the result is ( = 0.20829804e + 0.006858e 2 . In the case of 
random field disorder (short-ranged force-force correlations) ( — |. Both results compare 
well with numerical simulations. 

One should also note that (ITSl) contains a rather peculiar "anomalous term", namely 
R"'(0 + ) 2 , which only appears after the occurrence of the cusp. This term is in general hard 
to get, since the calculation naturally gives factors of -R'"(0), which are by parity, and not 
R"'(0 + ) = —R"'(0~). Several procedures to overcome these apparent ambiguities have been 
developed ll26ll . Supersymmetry will allow for another prescription, as will be discussed 
below. 

3. Supersymmetry and its Breaking 

Another way to average over disorder is to use additional fermionic degrees of freedom. It 
is more commonly referred to as the supersymmetric method. Supersymmetry is manifest 
using one copy of the system, where it immediately leads to dimensional reduction, as we 
show below. However it can not account for the non-trivial physics due to the appearance 
of the cusp and the corresponding breakdown of dimensional reduction, and supersymmetry. 
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This is possible when considering ti ^ 1 copies of the system, with n = 2 being completely 
legitimate. Here we give a general formulation, in which one can either discard the fermionic 
degrees of freedom, thus reconstructing a replica formulation at n = 0, or set n to e.g. n = 1, 
thus exploring supersymmetry. 
Define 

n f 1 2 

'ri[u a da,V} = ^ J - (Vu a (x)f + V(x,u a (x)) + ja{x)u a (x) . (16) 



a=l 



Then the normalized generating function of correlation functions for a given disorder V is 

jYla V l U a}e-^ HlUa ' ja ' V] 



(17) 



In the limit of T — > only configurations which minimize the energy survive; these 
configurations satisfy <5 ^"'^'^ = 0, of which we want to insert a ^-distribution in the path- 
integral. This has to be accompanied by a factor of det ^| , such that the integral 
over this configuration is normalized to 1, and supposing only a single configuration, the 
denominator can be dropped, leading to 



m = J \{v[u a }5 



5H[u a J ai V] 



det 



5 2 H[u a ,0,V] 
8u a (x)5u a (y) 



(18) 



5u a (x) 

Note that we neglect problems due to multiple minima, maxima, or saddle points. These 
configurations are incorrectly contained in (fTSl) . and are usually blamed for the failure of the 
supersymmetry approach. We will comment on this point later. For the moment, we continue 
with (fTSl and see how far we can get. Using an imaginary auxiliary field u(x) and two 
anticommuting Grassmann fields ijj(x) and ip(x) (per replica), this can be written as 

]JV[u a ]V[u a ]V$ a ]V[ip a ] 

5H[u a ,j a ,V] , - . .8 2 H[u a ,j a ,V] 



cxp 



- / u n (x 



5u a (x) r " v ' 5u a (x)5u a (y) 
Averaging over disorder yields with the force-force correlator A(u) : 



(19) 



-R"(u) 



S[u a ,U a ,4j a ,tfj a ,j a ] 



E 



J^J^>[w ]^[w a ]^[V'a]^ > [V'o] ex P {~<S[ u a, U a , tjj a , 1p a , j a ) ) 
U a (x){-V 2 U a {x) +J a {x)) +4> a {x){-V 2 )lp a {x) 

^ / ^tt a (a;)A(M a (x) - u h (x))u b {x) 

a,b ^ x 

+ -^a(x)lp a (x)A"(u a (x) - U b (x))4> b {x)lp b {x) 



(20) 



- u a (x)A'(u a (x) - u b (x))ip b (x)ilj b (x) 

We first analyze n = 1. Suppose that A(w) is even and analytic to start with, then only the 
following terms survive from (l20l) 



S Snsj [u } u } ^,ij,j} = / u(x)(-V 2 u(x)+3(x))+^(x)(-V 2 )^(x) - -u(x)A(0)u(x) (21) 
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(We have used that tp 2 = ip 2 = to get rid of the 4-fermion-term.) This action possesses a 
supersymmetry, which is most manifest when grouping terms together into a superfield 

u(x, e, e) = u(x) + &ip(x) + ij>(x)G + em(x) . (22) 

The action (f2~Tb can then be written with the SuperLaplacian A s as 

Ssusy = [ dede / u(x, e, e)(A 8 )u(x, e, e) , a, : = v 2 



and is invariant under the action of the supergenerators 



Q := 



■ev 



Q := x-=- + 



-ev . 



(23) 



(24) 



d& A(0) ' ^ " 96 A(0) 

Since "bosons" u and u, and "fermions" -0 and ip only appear to quadratic order, all 
expectation values are trivially Gaussian. Especially is 

A(0) 



UkU-k 



(k 



2\2 



(25) 



which is the result cited in ©, recalling that A(u) = —R"(u). Thus [u(x) — u{y)\ ~ 
A(0)|x — y\ 4 ^ d , which should be compared to the thermal average ([u(x) — u(y)} 2 ) ~ 
T\x — y\ 2 ~ d . Since both theories are Gaussian, the only difference is an appearant shift in 
the dimension of the system from d to d — 2. This is usually referred to as dimensional 
reduction. 

For more than n = 1 replicas, the theory is richer, and we will recover the 
renormalization of A(u) itself. To this purpose, write 



S[il a ,U a} ljj a} l/j a ,j a ] = ^ / \Ua(x)(-V 2 U a (x) + J a (x)) + lfj a (x)(-V 2 )lpa{ 

J x 



X) 



E 



-u a (x)A(0)u a (x) 



-u a (x)A(u a (x) - u b (x))u b {x) 

+ ^i> a {x)lpa{x)A"(u a (x) - U b (x))tfj b (x)lfj b (x) 

- u a (x)A'(u a (x) - u b (x))i/j b (x)i[) b (x) 



(26) 



Corrections to A(u) are easily constructed by remarking that the interaction term quadratic in 
u is almost identical to the treatment of the dynamics in the static limit (i.e. after integration 
over times) 



V__ -V 



+ 





(27) 



where an arrow indicates the correlation-function, a 
leads to (in the order given above) 



— y = (u(x)u(y)) = C(x - y). This 
8A(u) = [-A(u)A"(u) - A\u) 2 + A"(«)A(0)] / C(x - y f (28) 

J x—y 
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where the last term (being odd in u) vanishes. Note that this reproduces the non-linear terms 
inO. 

A non-trivial ingredient is the cancellation of the acausal loop in the dynamics (the 
"sloop", or 3-replica term in the replica formulation) ll26ll . This is provided by taking two 
terms proportional to u a A'(u a — u b )tjj b tp b , and contracting all fermions: 

^OvD =o - <29) 

since the fermionic loop (oriented wiggly line in the second diagram) contributes a factor of 
-1. 

One can treat the interacting theory completely in a superspace formulation. The action 

is 

a Je,eJx 

-\Y,I I- I- R(u a (x,e,e)-u b (x,e',e')). qo) 

2 a+b Jx Je,e J&,& 

Thus non-locality in replica-space or in time is replaced by non-locality in superspace, 
or more precisely in its anticommuting component. Corrections to R(u) all stem from 
"superdiagrams", which result into bilobal interactions in superspace, not trilocal, or higher. 
The latter find their equivalent in 3-local terms in replica-space in the replica-formulation, and 
3-local terms in time, in the dynamic formulation. 

Supersymmetry is broken, once A(0) changes. However, a new, shall we call it "effective 
supersymmetry", or "scale-dependent supersymmetry" appears, in which the parameter A (0), 
which appears in the Susy-transformation, changes with scale, according to equation (fl4l) . 

An interesting question is, whether anomalous terms, proportional to R"'(0 + ) can be 
recovered from the supersymmetric formulation. We show now, that this can indeed be done, 
in a very elegant way. The trick is to shift the disorder V(u) which appears in <f]"9t for the 
bosonic part, by a small amount S for the fermionic part, and to take the limit of 5 — ► at the 
end. This modifies d20b to 

= J2 Mx)(-V 2 u a {x)+j a {x))+t(j a {x){-V 2 )Mx) 

+ ^ / \]-Ua{x)A(u a (x) - U b (x))u b (x) + ^ a (x)ll) a (x)A"(u a (x) - U b (x))i) b (x)i) b (x) 



0,6 



(3D 



- u a (x)A'(8 + u a (x) - u b (x))ip b (x)i! b (x) 

I [u a (x)(-V 2 U a (x) + j a ( X )) + Va(x)(-V 2 )^ a (x) + ^ a (x) 2 A(0) 

- u a (x)A'(5)ip a (x)ip a (x) 

+ X/ / \^u a {x)A(u a (x) - u b (x))u b (x) + )-^ a (x)4) a (x)A''{u a (x) - u b (x))ifj b (x)ifj b (x) 

a^b -* x 
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- u a (x)lS!(5 + u a (x) - u b (x))4) b (x)i) h (x) . (32) 

Now the term with a = b is also well defined; the last term of (OTT) . as made explicit in (l3~2"l) . 
is of the form: 

/ ^2 -^a(x)A'(5)ip a (x)ilj a (x) - ^2 u a (x)A'(5 + u a (x) - u h {x))^) h {x)^ b {x) . (33) 

To demonstrate how this can be put into use, let us calculate the contribution to the 2-point 
function at 1-loop order, which is naively ambiguous lfT4irn?ll2"6l : 

<W P (u a (q)u a (-q)) = • ■ » --■ 





(q 2 + m 2 ) 2 A ^ J p (p + q/2) 2 + m 2 (p - q/2) 2 + m 2 ' (34) 
Note that the minus-sign comes from the closed fermion loop. In the limit of 5 — > 0, this gives 

5n oop M?KN) = - (?2+ 1 m2)2 A '(° + ) 2 S p {p + qm 2 + m 2 {p-q/2) 2 + m 2 ' ° 5) 

Also note that this would equivalently work for a iV-component field u — {u a } ,a = 1 . . . N, a 
case which poses additional difficulties, since derivatives have to be taken in a given direction 

Even higher correlation functions are immediate. E.g. is 

A^(0+) 4 /• 1 

1 '* ^ 64 ^(g2 +m 2)4 ■ ^ 

This has first been obtained, using straightforward dynamical perturbation theory in ll24ll . 
While the principle is simple, the calculations are actually very cumbersome. The sloop- 
method [26 1 is another efficient approach. (See also PUll .) 

Let us finally mention that the method correctly calculates the anomalous term of the 
"Mercedes- star" diagram at 3-loop order 031 . 

5R(u) = \ [J2» 4 - B!"{ufE!"{^f] , (37) 

where the icon stands for the momentum-integral only. The correction to A(u) = —R"(u), 
and picking the term proportional to A"'(u) = —R^ (it) comes from 



5A(u) = ...+ I \ =[... + A'(0 + ) 2 A'(u)A'»] /J\ (38) 





Note that two of the vertices are at argument A'(0 + ). This can otherwise only be calculated 
within the sloop-method. 
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4. Relation between Supersymmetry Breaking, Functional RG, and Replica-Symmetry 
Breaking 

Another popular approach to disordered systems is the replica-variational method, invoking 
replica- symmetry breaking (RSB). This method has for the problem at hand been developed 
in P4l . It consists in making the replacement 

S ^u a (x)A(u a {x) - u b {x))u b {x) — > S ^u a (x)u b (x)a ab . (39) 

a,b a,b 

This approximation is valid in the limit of a field u with an infinite number N of components; 
so for the following discussion we have to restrain ourselves to that limit. The variational 
replica approach then makes an ansatz for a ab , with different correlations a ab between 
different pairs of replicas. Finally a variational scheme is used to find an optimal a ab . The 
result (in the case of long-range correlated disorder, where the comparison can be made 123) 
is a hierarchic matrix with an infinite number of different parameters, of the form 




The exact form, and how it can be parameterized by a continuous function [a] (z), < z < 1 
is not of importance for the following. 

What is important is that RSB appears exactly at the same moment as in the functional 
RG the cusp appears ll35lfT8ll271l . Moreover, there is a precise relation between the 2-point 
function calculated by the RSB and FRG methods |[T%1I271 . The latter do not coincide, since 
the calculations are implicitly done in zero external field (RSB) and vanishing external field 
(FRG), leading to physically different situations (as for a standard ferro-magnet). 

In section|3]we have shown, that this is also the moment, when supersymmetry is broken. 
While the treatment there was for a 1 -component field, the conclusion is the same for an N- 
component field, and persists in the limit of iV — > oo. We can thus conclude that the breaking 
of supersymmetry, of replica symmetry, and the appearance of the cusp are all but different 
manifestations of the same underlying physical principle: the appearance of multiple minima. 
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